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An analysis of the vector meson spectrum from lattice QCD 

W. Armour * 

Department of Physics, University of Wales Swansea, Swansea, Wales, SA2 8PP 

We re-analyse meson sector data from the CP-PACS collaboration's dynamical simulations [XJ. Our analysis 
uses several different approaches, and compares the standard naive linear fit with the Adelaide Anzatz. We find 
that setting the scale using the J parameter gives remarkable agreement among data sets. Our predictions for 
the p and <j> masses have very small statistical errors, ~ 3 MeV, but the discrepancy between the different fitting 
approaches is ~ 40 MeV. 
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1. Introduction 

We study the chiral extrapolation of the vec- 
tor meson data from CP-PACS [I]. We do not 
have access to CP-PACS original data, so we pro- 
duced Gaussian distributions with central values 
and FWHM's equal to the quoted central values 
and errors respectively. Our data is uncorrelated 
throughout and we only fit to degenerate data, 

1. e. k;^ = K 2 val . The CP-PACS data used is 
from mean-field improved Wilson fermions with 
improved glue at four different (3 values. For each 
/? value there are four different « soa values, giv- 
ing sixteen independent ensembles. The physical 
volume was held fixed at La w 2.5 fm for the 
f3 = 1.80,1.95 and 2.10, but the (3 = 2.20 en- 
semble had a slightly smaller physical volume. A 
graphical overview of the CP-PACS data is given 
in Figure d We also fit to CP-PACS quenched 
data for comparison. 

2. Fitting analysis 

2.1. Summary of analysis techniques 

Our chiral extrapolation approach is based 
upon converting all masses into physical units 
prior to any extrapolation being performed. An 
alternative approach would be to extrapolate di- 
mensionless masses (in lattice units) PQ. Our 
method has the following two advantages: 

• Different ensemble's data can be combined 
together in a global fit. 



*The author wishes to thank PPARC for funding. 
C. Allton for his guidance, help and contributions, also 
the Adelaide group for there useful comments. 
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Figure 1. A Graphical overview of the CP-PACS 
data. Note that here the scale is set using the 
inverse lattice spacing determined from the J pa- 
rameter |2]. 

• Dimensionful mass predictions from lattice 
simulations are effectively mass ratios, and 
hence one would expect some of the system- 
atic and statistical errors to cancel. 

2.2. Fitting functions 

In our chiral extrapolations we use the naive 
linear fit, Eq.JJJ, as well as the Adelaide Anzatz, 
gj Eq.0. 



My(/3, K sca ; «W, «val) = 

Co + C 2 Mp S ({3, K soa ; K va i, K va i). 
My(/3, K soa ; K val , K val ) = 

Co + C2Mp S (/3, K soa ; «val, Kval) + 
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where the Self Energies are defined as 

£^ 7r (Mp S (/3, K soa ; K soa , Kvai)) and ££ w 



the 



^w( M Ps{0> K sca5 K soa , Kval))- AV(PS) 

vector(pseudo-scalar) meson mass. (3 and K soa 
refer to the sea structure (i.e. the gauge coupling 
and sea quark hopping parameter) and n va i refers 
to the valence quark hopping parameters. 

The Adelaide method relies on a parameter A ffW 
0. We use the value of A™ = 630 MeV taken 
from The value of A^ is highly constrained 
by the lightest data point in the My versus Mp S 
plot, and since the data used in jS] includes a 
much lighter point than in this study, we use its 
value of A^. 

2.3. Individual ensemble fits 

Our first method for obtaining physical masses 
involves determining fitting parameters, Co and 
C*2, for each of the 16 data sets and then perform- 
ing a continuum extrapolation to those. Figure[21 
displays the Co values and motivates a continuum 
extrapolation of the form in Eq.0. 
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Figure 2. An example of a continuum extrapo- 
lation performed to fitting parameters obtained 
from the sixteen different data sets. 



Co, 2 (a) — Co° 2 "' + A 0i 2 o,j- 

Table 1 lists the results for these fits. 



(3) 



2.4. Global fits 

Figure [3] plots the vector meson mass against 
the pseudo-scalar mass squared (in physical 
units) for all of the data from all of the ensembles. 
As can be seen the data lies on a near universal 
line. This motivates an analysis which combines 
all of the degenerate data into one global fit. 
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Figure 3. A plot of My versus Mj> s for the de- 
generate CP-PACS dataset showing that the data 
lies on a near universal line. Also plotted are the 
three experimental data points, and the unitary 
data from UKQCD gj. Note that the scale for 
the lattice data was set using the method of 0- 

To underline the near universal nature of the data 
we plot the relative spread of the data using EqEI 



spread = 



M 



data 



Co + C^-Mpg 



(4) 



Where Co + C 2 Mj> s is a simple fit through the 
data in Figure 3. In Figure 01 we see that the 
spread is at the 1% level or less. Note that this 
simple linear fit is not used in any further anal- 
ysis. Also note that the scale is set from J [2] 
and hence the data must pass through the (M#-* 
, M.x) point. When grouping the different en- 
sembles together we must allow for some (small) 
variation with the lattice spacing. We do this by 
modifying the linear and Adelaide fitting func- 
tions. By substituting Eq.© into Eqs.lP&EJ) 
respectively. Table 2 lists the results of these fits. 
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^jcont. 
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[GeV] 


[GeV/fm] 




[GeV- 1 ] 


[GeV-VH 




Linear-fit 


0.772 ± 3 


8±2 x 10" 2 


1.10/14 


0.473 ± 6 


-0.27 ±3 


1.86/14 


Adelaide- fit 


0.785 ± 3 


7±2 x 10" 2 


2.01/14 


0.462 ±6 


-0.26 ±3 


1.39/14 



Table 1: The coefficients obtained from the continuum extrapolation of the parameters obtained from 
the 16 individual fits using Eq.(3). 
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[GeV" 1 ] 


[GeV- 1 /fin] 




Linear-fit 


0.772 ± 3 


8±2 x 10" 2 


0.474 ±6 


-0.28 ±3 


43/76 


Adelaide-fit 


0.784 ± 3 


7.8 ± 16 x 10~ 2 


0.466 ±6 


-0.28 ±3 


58/76 



Table 2: The coefficients obtained from a global fit of all the My data against Mp s using the linear-O(a), 
Eq.(l), and Adelaide- O(a), Eq.(2), fits, both incorporating Eq.(5). 
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Figure 4. A plot showing the relative spread in 
My versus Mp s for the degenerate CP-PACS 
data set. 



3. Conclusions 

To conclude, Table [3] lists our mass estimates 
for the p and 4> mesons. We see that our inter- 
pretation of the Adelaide Anzatz underestimates 
M p presumably due to a poorly tuned A^ value. 
We also note the following: 

• Setting the scale using J gives remarkable 
agreement among data sets. 

• The (statistical) errors in the mass esti- 
mates are tiny. 

• The discrepancies between the various fit- 
ting procedures is much larger than the sta- 
tistical errors listed. 



• Note that for the global linear fit that incor- 
porates the Xa correction, we obtain an M p 
only 10 MeV above the experimental value. 

• The estimates of M p and M$ from this ap- 
proach are closer to the corresponding ex- 
perimental values than the quenched esti- 
mates. 





M p 


M 




[GeV] 


[GeV] 


Experiment 


0.770 


1.0194 


Quenched + Xo,2 


0.798 ±4 


0.988 ±5 


Global-Linear + Xo,2 


0.781 ±3 


0.995 ±2 


Global-Adelaide + Xq,2 


0.740 ± 3 





Table 3: Mass predictions. 
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